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1. Introduction 

Very-high-energy (VHE) gamma-ray sources are a fascinating probe into cosmic particle ac¬ 
celeration processes in the most extreme environments of the Milky Way Galaxy. There are a 
number of Galactic VHE source classes that have been detected including supernova remnants 
(SNRs), pulsar wind nebulae (PWNe), binary systems, and star forming regions, as well as yet 
unidentified sources of VHE emission. SNRs and PWNe, in particular, are the likely accelerators 
of cosmic rays of Galactic origins. PWNe may also be the source of the positron excess detected 
by Pamela[l]. 

In recent years the H.E.S.S. Galactic Plane Survey (HGPS) has been a fruitful probe of the 
central area of the Galaxy [2], detecting numerous new VHE sources, and probing diffuse Galactic 
VHE emission. The Cygnus region is a known region of VHE gamma ray emission, with emission 
detected by Milagro [3], HEGRA [4] [5], as well as VERITAS [6] [7] [8]. It also contains a number 
of potential gamma-ray emitters, identified by fhe EGRET [9] and Fermi-LAI GeV experimenfs 
[10] [11]. Mofivafed by fhe densify of pofenfial VHE gamma-ray emiffers, VERITAS undertook 
survey observations of fhe region from 2007-2009. These observafions, as well as follow-up ob¬ 
servations have been analyzed fo provide fhe deepesf look af fhe Cygnus region af TeV energies fo 
dale. 

2. VERITAS Survey and Follow-up Observations 

VERITAS is an imaging afmospheric Cherenkov telescope (lACT) array system locafed af fhe 
Ered Eawrence Whipple Observatory (EEWO) in soufhern Arizona (31°40' N, 110°57' W, 1.3km 
a.s.L). The array consisfs of four 12 m diameter telescopes wifh a field-of-view of 3.5°, which 
deled gamma rays af energies of 85 GeV fo > 30 TeV af an energy resolulion of 15-25%. The 
pointing accuracy of fhe array is belter lhan 50 arcsec and fhe angular resolution is < 0.1° wifh a 
68% conlainmenl radius af 1 TeV. Eour-lelescope observations began in 2007, and Ihere have been 
Ihree major upgrades in fhe subsequenl years. A lelescope move in 2009 improved fhe sensilivily 
of fhe array [12], and a Irigger upgrade in winter 2012 and a camera upgrade in summer 2012 have 
improved fhe sensilivily furlher and lowered fhe energy fhreshold of fhe array. 

Erom 2007 Ihi'ough 2009, VERITAS undertook a survey of a 15° x 5° area of fhe Cygnus 
region of fhe Galaxy, covering Galaclic longilude (1) from 67° to 82° and Galaclic lalilude (b) from 
— 1° fo 4°. The survey consisled of more lhan 140 hours of observations, reaching an average poinl- 
source sensilivily of belter lhan 4% Crab Nebula flux al energies above 200 GeV [13]. The base 
survey observations achieved a nearly uniform effective exposure of ~6-7 hours al every location 
over Ihe field of Ihe survey. Observations were scheduled to keep Ihe average zenilh angle around 
20° as much as possible in order to Iry to avoid Ihe higher energy Ihreshold associated wilh large 
zenilh angle observations. The full survey dalasel also includes over 150 hours of follow-up and 
poinl source observations in Ihe region (laken from 2008 Ihrough 2012). Eollow-up dala were laken 
using Ihe slandard VERITAS wobble observing slralegy for poinl sources [14]. 

The previous blind search analysis from Ihe survey plus follow-up observations yielded Iwo 
source detections and suggested possible VHE /-ray emission al several olher locations. There 
were Ihree regions observed in follow-up observations, shown in Eigure 1. Eour sources were 
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right ascension 



galactic longitude 



Figure 1: The radial acceptance corrected exposure map for the VERITAS observations, including follow¬ 
up observations. The scale on the color bar is in hours. The yellow box demarcates the core survey region, 
and the gray circles indicate the follow-up observation regions. 

detected: VER J2019-I-407 (G78.2-1-2.1, y-Cygni) is likely a SNR shock front interacting with am¬ 
bient material, and demonstrates hadronic emission [6], TeV J2032-I-413, the first unidentified TeV 
source, and likely PWN [8], and in fhe OBI region (coincidenf wifh MGRO J2019-I-368) VER 
J2016-I-371 (a poinf source associafed wifh CTB 87 and a likely PWN) and VER J2019-I-368 (a 
broad excess wifh some subsfrucfure suggesting the possibility of multiple sources) [7]. We will be 
able to update and expand on these results with the improvements in the analysis over the last few 
years. 

3. Analysis 

3.1 VERITAS Analysis 

The analysis of the Cygnus region data set begins in the same way as the standard VERITAS 
analysis with a calibration and cleaning of the shower images, followed by a principal component 
analysis of the shower images [15]. The trials factors of the analysis have been reduced by limiting 
the number of analyses taken with cuts. The analysis choices are motivated by the characteristics 
of know Galactic sources. The a priori cuts include a choice of two spectral cuts: one for mod¬ 
erate and one for hard sources. The spectral cuts take into account the telescope upgrades, and 
are optimized for the spectral index and strength of the potential source, and for different energy 
thresholds. There is also an analysis with known sources are defined according fo fhe published 
VERITAS analysis values for exfension. The backgrounds are esfimafed using fhe ring background 
mefhod applied af each sky map bin following fhe mefhods of [16]. 

In order fo opfimize fhe analysis sensifivify while limifing fhe number of frials over fhe large 
skymap we underlook an analylical calculalion following [17] fo defermine fhe optimum bin size 
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for the analysis. Because Galactic sources are typically extended, we assumed a two-dimensional 
Gaussian with a la (1 standard deviation) containment diameter of 0.5° on a 10% Crab source with 
15 hours of exposure. The optimum bin size is linearly dependent on the width of the Gaussian and 
inversely proportional to the cube root of the intensity. The optimum bin size was found to be just 
a bit under 0.1°; thus, in this analysis the results are re-binned to 0.1° prior to further analysis. 

Furthermore, in order to fully account for the trials in the analysis we are performing a Monte 
Carlo simulation [18]. In each step of the simulation a map of randomly generated background 
gamma-rays is populated according to the dead-time and radial acceptance corrected exposure 
map. Then a ring background analysis is done on each bin of the map, and the highest significance 
bin and a random bin significance are recorded. In this analysis the random bin is picked to be in 
the core survey region, and not in the edges of the sky map. This is repeated about 10^ times to 
get large the statistics to fit to high significances. The frials facfor of fhe VERITAS analysis will be 
faken info accounf in fhe final resulfs. 

3.2 Fermi-hAIT Analysis 



Figure 2: A Fermi counts map of the Cygnus (>1 GeV) region using nearly six years of data, Fermi 3FGL 
sources are labeled. The square denotes the VERITAS sky survey region. The color bar scale is in square 
root of counts. 


An analysis of six years of Fermi-LPTV data (Pass 7) in the Cygnus region in the energy range 
from 1 - 300 GeV has been undertaken using the publicly available EAT data tools ^ and the binned 
analysis technique. Since this is a large region the analysis was initially broken up into three regions 
of 10° radii centered at (l,b) of (79.5°, 1.5°), (74.5°, 1.5°), and (69.5,1.5°). Sources outside of 7° 
and with low significances in fhe 3EGE catalog (<5a) [10] were initially fixed. The resulfs of fhese 
fifs were fhen combined fo a fofal fif of a 15° region cenfered af fhe middle of fhe region. In each 
of fhese fifs fhe exfended sources in fhe region were held fixed and modeled using fhe femplafes 
provided by fhe Fenni-LKV collaborafion. A model and a counfs map (Eigure 2) were compared 

*7erra!-LAT tools v9r33pO. They are available at http://fermi.gsfc.nasa.gov/ssc/data/analysis/ 
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to check for any residuals, and the region is well fit by the model. The results from the fit are 
summarized in Table 1 . 

There are 25 3FGL catalog sources in the VERITAS survey Cygnus region (labeled in Fig¬ 
ure 2) including two that are low significance (< 5a) in fhe cafalog, and fwo fhaf are exfended 
sources (Gamma Cygni, and fhe Cygnus Cocoon). The 21 remaining sources were fil using 
fhe binned analysis described above. Two of fhem are of pofenfial exfra-galacfic origin: 3FGL 
J2015.6-t3709 (associafed wifh MG2 J201534-t3710, an FSQR) and 3FGF J2018.5-t3851 (asso- 
ciafed wifh TXS 2016-1-386, an acfive galaxy of uncerfain fype). Of fhe sources wifh Galacfic 
associafions: 3FGF J2032.2-I-4126, is associated wifh fhe FAT pulsar J2032-I-4127, and is fenfa- 
fively associafed wifh fhe unidentified TeV source J2032-I-4130, which could have a PWN asso¬ 
ciation, and 3FGF J2021.1-1-3651 and 3FGF J2017.9-1-3627 are associated wifh fhe region of fhe 
extended Milagro source MGRO J2019-I-37. Two of fhe remaining sources (3FGF J1952.9-I-3253 
and 3FGF J2021.5-1-4026) are associafed wifh pulsars, and fhe remaining 14 are unidentified. All of 
fhe sources fhaf were fil in fhe analysis of fhe Fermi-LPH dafa were sfudied for promising charac- 
ferisfics for pofenfial TeV emission such as a sfrong GeV defecfion significance, a hard speclrum, 
and a high exlrapolaled flux. We will cross correlafe fhe resulfs and maps of fhis analysis wifh fhe 
VERITAS resulfs. 

Furthermore, Ihere are nine sources wilhin fhe VERITAS survey region in fhe Firsf Fermi- 
FAT Cafalog of Sources above 10 GeV (designated IFHF) [11], which uses dafa from fhe firsf 
Ihree years of fhe Fermi mission. Included among fhese sources are fwo sources associafed wifh 
acfive galacfic nuclei, four pulsars wifh pulsafion above 10 GeV (fwo of which are associafed wifh 
fhe TeV sources MGRO J2019-I-37 and TeV J2032-I-4130), one GeV source associated wifh fhe 
SNR G078.2-I-02.1, one associated wifh fhe slar forming region of fhe Cygnus Cocoon, and one 
unidenlified source. Due lo fhe longer observafion window of our analysis and fhe longer time 
window of fhe 3FGF cafalog, 3FGF cafalog sources were used lo populale fhe binned analysis 
inpuls. 


4. Conclusions and Outlook 

The VERITAS survey of fhe Cygnus region serves lo enhance fhe knowledge aboul fhe popu¬ 
lations of VHE gamma-ray emillers, and pofenfial VHE gamma-ray emitters in our Galaxy, using 
over 300 hours of observafions in fhe region from 67° lo 82° Galactic longilude and from —1° lo 
4° Galacfic lalilude. This region has been defermined lo be acfive in VHE gamma rays wifh four 
sources detected by VERITAS and is furthermore a region for active inquiry at GeV energies with 
a large number of unidentified sources. The new VERITAS analysis in the region has increased 
sensitivity, and can help untangle the origins of gamma-ray emission in the region. 

In the future, there will likely be further VERITAS observations in the region as a result of this 
analysis. The cross correlation of the Fermi and VERITAS sky maps will provide strong motivation 
for low significance source candidales, and will allow for fhe possible idenlificafion or association 
of unidenlified GeV sources in fhe region. Furlhermore HAWC will have good sensilivify in fhe 
fhe region and VERITAS will work fo follow-up on pofenfial new sources delecled wifh HAWC. 
Togelher fhis will enable a grealer underslanding of fhe region al fhe highesl energies. 
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3FGL name 

Association(s) 

TS 

Flux (1 - 300 GeV) 

Spectral shape 

Extrapolated Flux 
(300 GeV -1 TeV) 

(cra-Vl) 

J1952.9+3253 

PSR J1952+3252 

IFHL J1953.3+3251 

5973 

L96E-08+/-4.17E-10 

broken power law 

L19E-09 

J2004.4+3338 

IFHL J2004.4+3339 

246 

3.48E-09+/-3.04E-10 

power law 

6.97E-13 

J2015.6+3709 

MG2 J201534+3710 

IFHL J2015.8+3710 

908 

9.06E-09+/-4.30E-10 

log parabola 

LllE-13 

J2017.9+3627 

MGRO J2019+37 

667 

9.84E-09+/-5.16E-10 

log parabola 

L72E-15 

J2018.5+3851 

TXS 2016+386 

IFHL J2018.3+3851 

76 

L88E-09+/-3.09E-10 

power law 

7.39E-13 

J2018.6+4213 

PSR J2021+3651 

34 

L41E-09+/-2.90E-10 

power law 

L27E-12 

J2021.1+3651 

MGRO J2019+37 

IFHL J202L0+3651 

23561 

6.84E-08+/-8.04E-10 

broken power law 

4.15E-10 

J2021.5+4026 

LAT PSR J2021+4026 

46339 

L13E-07+/-9.45E-10 

broken power law 

2.79E-10 

J2022.2+3840 


98 

3.35E-09+/-3.90E-10 

log parabola 

4.50E-20 

J2023.5+4126 


84 

3.29E-09+/-4.25E-10 

log parabola 

3.52E-14 

J2024.6+3747 


28 

L37E-09+/-4.78E-10 

log parabola 

3.69E-15 

J2026.8+4003 

LAT PSR J2032+4127 

68 

3.19E-09+/-4.46E-10 

log parabola 

L13E-14 

J2032.2+4126 

TeV J2032+4130 

4717 

2.45E-08+/-5.74E-10 

broken power law 

L23E-09 


IFHL J2032.1+4125 





J2032.5+3921 


132 

4.38E-09+/-4.32E-10 

log parabola 

6.70E-16 

J2032.5+4032 


116 

4.20E-09+/-6.33E-10 

log parabola 

7.13E-22 

J2034.6+4302 


69 

2.65E-09+/-5.01E-10 

log parabola 

L60E-18 

J2036.8+4234C 


75 

2.52E-09+/-4.38E-10 

power law 

4.53E-14 

J2037.4+4132C 


62 

2.42E-09+/-4.24E-10 

power law 

L19E-13 

J2038.4+4212 


129 

3.64E-09+/-4.68E-10 

log parabola 

3.20E-14 

J2039.4+4111 


75 

2.71E-09+/-3.98E-10 

log parabola 

L97E-14 

J2042.4+4209 


76 

2.68E-09+/-3.82E-10 

log parabola 

4.70E-14 


Table 1: The table summarizes the results of the fit of six years of Fermi-LAT data in the Cygnus region. 
The 21 sources used for the binned analysis are listed together with its associated source(s) if applicable 
and a test statistic (TS) value, the flux from 1 - 300 GeV, and an extrapolated flux at 300 GeV - 1 TeV is 
calculated using the spectral shape noted in the table. 
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